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Abstract 

A two-dimensional generalized oscillator with time-dependent parameters is considered 
to study the two-mode squeezing phenomena. Specific choices of the parameters are used 
to determine the dispersion matrix and analytic expressions, in terms of standard hermite 
polynomials, of the wavefunctions and photon distributions. 


1 Introduction 

In the middle of the sixties and beginning of the seventies a set of quantum states of the elec- 
tromagnetic field were observed which have less uncertainty in one quadrature than a coherent 
state [1-3]. These one-mode squeezed states have generated big expectations in optical commu- 
nication systems [4]. In some quantized fields, the interaction hamiltonians occur only between 
pairs of modes and then to understand the main features of the system, one restricts to study 
one and two noimal modes . In the last decade two- mode squeezing phenomena have attracted 
attention to study properties of noise and correlations [5-8]. Recently the accidental degeneracy 
of a two-dimensional (2-D) harmonic oscillator with frequency plus an interaction proportional 
to the z-th projection of the angular momentum was studied [9]. This system was called the 
generalized 2-D harmonic oscillator because presents a bigger accidental degeneracy depending on 
the strength A of the angular momentum interaction. This model was generalized [10] to include 
time-dependent parameters, m = m 0 f{t) and A = u; 0 A 0 (f). If we take /( 0) = A 0 (0) = 1; the 
hamiltonian, for t = 0, represents a charged particle moving in a constant magnetic field . 

The aim of this work is to study two-mode squeezing phenomena with this model because it 
demonstrates the change of dispersions due to variation of the mass and coupling constant during 
the evolution. In the framework of quantum optics the hamiltonian is built by: the operator 
(1 / f T that causes a time-dependent, exchange of kinetic and potential energies within 

each mode, the interaction 1/2(1// — f)((P • $ + a • a), which describes a degenerate two-photon 
interaction; and the potential «Ao(t)(a 2 a i — a{ 02 ) 1 that is a mode mixing operator. 

The solution of the corresponding time dependent Schroedinger equation is obtained through 
the theory of integrals of motion [11]. By means of Noether’s theorem, using a special variation 
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[10], we construct the linear time dependent integrals of the motion. The resulting quantum 
invariants are given in terms of the positions and momenta operators [10,11] by 

P(t) = A xp+ A 29, Q(t) = A 3 p+ A 4 g, ( 1 ) 

with the initial conditions P{ 0) = p and Q(0) = q, so that the 2 x 2_mat rices previously introduced 
satisfy Aj(0) = A4(0) = / 2 and A3(0) = A2(0) = 0. The operators A(t) = 1/ y/2 [Q(t)/l + il/hP(t)] 
and its hermitean conjugate, can be constructed with the matrices 

A p = yA3 + — Ai, A g = yA 4 + — A 2 , ( 2 ) 


with l = defining the oscillator length. These integrals of motion also are given in terms 

of the creation and annihilation photon operators 

A(t) = Mi a + M 2 a*, A\t) = M$a + . (3) 


With the initial conditions A( 0) = a and A f (0) = a*, the matrices defined in (3) comply with 
Mj(0) = M 4 (0) = h and M 3 (0) = M 2 (0) = 0. The A fc ’s, M fc ’s, A p and A, are entries of symplectic 
matrices in four dimensions because the invariants (1) and (3) satisfy the commutation relations 
of Heisenberg- Weyl algebras. 

In the present work, we study the behavior of the model for Ao(t) an arbitrary function of time 
and considering two kinds of varying masses, i.e., two choices for the function f(t), namely: 


f(t) = exp( 7 1 ) ; 

f(t) = J cosh 2 Clot 

l {Cl 0 (t — TjsinhfloT 1 + coshCloT} 2 

For these two cases the A* matrices take the general form 


( cos 0 sin# A 

x k = p k n- p k y_ siri6 COS 0j; 



( 4 ) 

t < 0 

0 <t<T . 
T < t 

( 5 ) 


= 1, 2, 3, 4; (6) 


where the definition 9 = Jq lli 0 X 0 (t)(It was used. The analytic expressions for the p k ’s functions 
are given in Ref. [10]. In the next sections we determine the coherent and Fock-like states, the 
photon distributions and the dispersion matrices in terms of these /ujt’s. 


2 Squeezed Coherent and Fock States 

The coherent-like states are obtained by solving the differential equation A(t)$ Q {q,t) = 0 with 
A(t) given in Eq. (3). This solution yields the vacuum state of the physical system, and its phase 
is chosen to guarantee that satisfies the time dependent Schroedinger equation. The expression 
for the ground state wavefunction is 

1 f * Pq — -.1 

7 SF^ exp i^‘ 9 } ‘ 


$0 (q,t) = 


( 7 ) 



To get the last expression the relation (2) was used and the functions fi p = and 

Vq — 2 (jf/^ + 1 ^ 4 ) were defined. To obtain the general expression for the eigenstates in the 

coordinate representation one needs to apply the unitary operator D(a) = exp{<5 • A* - a* • A], 
which is an invariant, to the vacuum wavefunction (7), i.e., 

*„(?,«)= exp {-ld. + I^ 5 . 5+ _i. ? -R 5 } w>() . (8) 

These are expressed in terms of multi-dimensional Hermite polynomials [12] through the relation 

( 1 _ \ 00 A '*" 1 / V * n » /„ \ 

<r + t,Q-R 7 )= Y, -fr , HS(-Rf) • (9) 

/ rn , n 2 =0 U 1‘ '*2- ' w t 


Substituting the last expression into (8) and using the form of the coherent-like states in the 
Fock-like representation, we get the Fock-like eigenstates in the coordinate representation: 


(q\nxn 2 ) = $ 0 (q,t) H 



711,712 



( 10 ) 


These multi-dimensional Hermite polynomials are rewritten as a product of two standard one- 
dimensional Hermite polynomials [12] as follows: 


„ I f / i \ / u*\(ni+n2)/2 / I \ 

H ”‘~ (~h^ r R V = Vifr) H "‘ (vi%^ lcos '' qi + sin0 *0 

* (ti4u i_ sinS q ‘ + cos6 ,2 0 ’ (u) 

where we use the explicit expression of matrix R. These Fock (10) and coherent (8) -like states 
represent squeezed and correlated eigenstates of the system as it will be shown further. 


3 Propagator 


The propagator in the coherent state representation is given by the matrix elements of the evolution 
operator U(t), which will be obtained by means of the theory of time dependent integrals of 
motion [11]. If I(t) is an integral of motion then satisfies I(t)U(t) — U(t)I( 0). Taking its matrix 
elements with respect to the coherent states, we get a linear system of differential equations, which 
can be solved. Thus the propagator takes the form 


G(<r,7,t) = exp ( l7|2/2) exp + S-M^+ . 

For the cases (4) and (5), the following relations are satisfied 

\JdetMx = ~ —^pj 9 i> ! M 2 = "^2^ (^Vq + = 92^2 , 


( 12 ) 


(13) 
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( 14 ) 


Substituting these relations into the Eq.(12) we get the propagator, which through Eq. (9) can 
be expressed in terms of multi-dimensional Hermite polynomials. If we compare with the power 
series expansion of the propagator we get the probability amplitude for having n x and n 2 photons 
in the coherent-like state |7, t), t.e., 

( n ." a |7,t) = Wt ),17^ ex P Hf- + ' (15) 

By means of the Eq.(12) this amplitude can be rewritten in terms of standard Hermite poly- 
nomials [12]. The squared absolute value of this amplitude yields the photon distribution function 
of the system, W„ inj (7,«) = |(nin 2 |7,t)| 2 . This will let us calculate, at least formally, the mean, 
(N k ), and the mean squared fluctuation of the number of photons, (A N k ) 2 , in direction fc, which 
are present in the coherent state |7, t). The expectation values of N k and N k are evaluated directly 
using the expressions of the creation and annihilation photon operators in terms of the integrals 
of the motion (3), and the commutation properties for these invariants. For the vacuum state one 
has 2 

(N k ) = ^{(^1 “ l 1 *) 2 + (jn Q u 0 fi 3 -I- m Pi) } > ( 16 ) 

{N%) = - Ha) 2 + + -777M2) 2 } 

{(«+*) 2 + ("W<3 ' (17) 

With these expressions, we evaluate the ratio of the mean squared fluctuation (A N k ) 2 and the 
mean number of photons ( N k ), which determines the nature of the distribution function of the 
system: 

= 5 { (,i| + " 4)2 + ( m ° u '°'‘ 3 “ 1 <18) 

For the cases (4) and (5) the ratio is greater than one when t > 0, which implies that we have 
a super- Poissonian photon distribution function. For t = 0, there is a discontinuity in the ratio, 
which is obtained by comparing the following limiting procedures: making t — ► 0 and then a — ► 0, 
and conversely. 

4 Dispersion Matrices 

The dispersion matrix can be written in terms of 2 x 2 matrices characterizing the dispersions 
in the positions and momenta operators and the correlation between them. Besides for the cases 
under study, due to (6), they take the form 

= h • (19) 
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^ (<) = " 2 + 12 ’ 

° q q (t) = 5— 77(^1 + (mowto) 2 /^) I 2 . 

The corresponding correlation matrices for the creation and annihilation operators 
immediately from the last expressions; they are given by 


( 20 ) 

( 21 ) 

are obtained 


= j{" ? " A + (ra »“» )2 ^ " - 2i ( 


1 


(m 0 w 0 ) 2 


\rn 0 u > 0 


H1H2 + m 0 u>on 3 fj, 4 


)} 


( 22 ) 


+ , ( 23 ) 

= ~ "* + (m " u - ,o)2 '’ 2 _ (m ^ o) 2 + 2i (^;'‘ 1 '' 2 + } • ( 24 ) 





(•) (b) 

Fig. 1. Dispersion and correlation matrices behavior in positions and momenta 
space for the studied cases in this paper: (a) corresponds to Eq.(4), and (b), to Eq.(5). 
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The behavior of the dispersion matrices is illustrated in Fig. 1. For the case (4), we choose the 
parameters 7 = 0.1 and m 0 = w 0 = 1. It is seen that there is squeezing for the coordinates and 
stretching for the momenta. Also one notes that is a negative function and therefore there 
are one-mode correlations between the coordinates and the momenta. If we reverse the sign of 7, 
the roles between the dispersion for coordinates and momenta are interchanged, and a M becomes 
positive. In the case (5), we use the parameters flo == 0.15, T = 10, and mo = u>o — 1* In spite 
of the mass is different that in the previous example, the general trends are similar. For example, 
the <jpp is an increasing function of time starting from its minimum value at t < 0, and there is 
squeezing for the o qq . The main difference appears in the correlation in this case, it can be 
positive for large times, while in the previous one is negative or zero for any time. 
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